Lamellar order, microphase structures and glassy phase in a field theoretic model for 

charged colloids 



Marco Tarzia a and Antonio Coniglio a ' b 
a Dipartimento di Scienze Fisiche and INFN sezione di Napoli, Universitd degli Studi di Napoli "Federico IF', 
Complesso Universitario di Monte Sant'Angelo, via Cinthia, 80126 Napoli, Italy and 

b Coherentia CNR-INFM 
(Dated: February 6, 2008) 

In this paper we present a detailed analytical study of the phase diagram and of the structural 
properties of a field theoretic model with a short-range attraction and a competing long-range 
screened repulsion. We provide a full derivation and expanded discussion and digression on results 
previously reported briefly in M. Tarzia and A. Coniglio, Phys. Rev. Lett. 96, 075702 (2006). The 
model contains the essential features of the effective interaction potential among charged colloids in 
polymeric solutions. We employ the self-consistent Hartree approximation and a replica approach, 
and we show that varying the parameters of the repulsive potential and the temperature yields a 
phase coexistence, a lamellar and a glassy phase. Our results suggest that the cluster phase observed 
in charged colloids might be the signature of an underlying equilibrium lamellar phase, hidden on 
experimental time scales, and emphasize that the formation of microphase structures may play a 
prominent role in the process of colloidal gelation. 
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I. INTRODUCTION 

Colloidal suspensions are solutions of solid (or liquid) 
mesoscopic particles immersed into another substance [l[ . 
These systems, like blood, proteins in water, milk, inks, 
or paints, are ubiquitous in our everyday life and are 
extremely important in biology and industry. Due to 
their potential applications for designing new materi- 
als with a wide range of viscoelastic properties, in the 
last few years there has been much interest in the role 
of the inter-particle potential on controlling the struc- 
tural and dynamical properties of colloidal systems. By 
appropriately varying some control parameters (such as 
the composition of the solvent, the coating of the parti- 
cles, the concentration of the polymers into the solvent, 
. . . ) the effective interaction potential between colloids 
can be suitably tuned in the experiments. It is possi- 
ble to realize a hard-sphere system jsj; by adding non 
adsorbing polymers the hard-sphere interaction can be 
complemented by a short-range attraction, induced by 
depletion forces |3J . Recent experimental works outlined 
that in some cases a residual net charge on the surface 
of colloidal particles may be present [j, Q, thereby in- 
ducing a long-range electrostatic repulsion screened by 
the presence of ions in the solution. The resulting ef- 
fective interaction is therefore given by a hard-core term 
accounting for the excluded volume, a depletion-induced 
narrow attractive shell and a long-range repulsive shoul- 
der. This kind of potential is well approximated by the 
DLVO (after Derjaguin, Landau, Verwey and Overbeek) 
potential Q . Interestingly enough, the parameters of the 
potential can be suitably tuned in the experiments. The 
depth of the attractive shell is controlled by the concen- 
tration of the polymers, <p p , which therefore plays the 
role of an inverse temperature, whereas its range is pro- 
portional to the ratio between the radius of the colloids, 



a, and the gyration radius of the polymers, R g (typically 
one has 0.05 < R g /a < 0.2). On the other hand, the ad- 
dition of salt increases the number of ions in the solution 
responsible for the screening of the electrostatic repul- 
sion, thereby reducing the amplitude and increasing the 
screening length of the Yukawa potential. 

Charged colloids in polymeric solutions have recently 
raised a lot of interest .both from an ex per imental 0, 0] 
and a theoretical @, \& \M, EH, E3, 0, 13 point of view. 
In these systems the competition between attractive and 
repulsive interactions on different length scales stabilizes 
the formation of aggregates of an optimal size and shape 
(cluster phase), characterized by a peak of the struc- 
ture factor around a typical wave vector, k m . Experi- 
mentally, such cluster phase, made up by approximately 
monodisperse equilibrium aggregates, can be clearly ob- 
served using confocal microscopy at low volume fraction 
and low temperature (high attraction strength) 0, H, 0] . 
By appropriately tuning the control parameters (i.e., in- 
creasing the volume fraction or decreasing the temper- 
ature) the system progressively evolves toward a gel- 
like non-ergodic disordered state 0, 0] (colloidal gela- 
tion), where structural arrest occurs. Although intensely 
studied both experimentally and numerically, a theoret- 
ical understanding these phenomena is still lacking: the 
mechanisms inducing colloidal gelation are still unknown 
and many gaps remain in our present knowledge of the 
equilibrium phase diagram of these systems. InRef. [H, it 
has been proposed that colloidal gelation is related to the 
formation of a Wigner glass, whose blocks are compact 
and thermodynamically stable clusters, with a residual 
long-range repulsion between them. In Ref. [lOj it was 
instead highlighted the percolative nature of the struc- 
tural arrest in charged colloids, suggesting that perco- 
lation of a spanning network of clusters with long living 
bonds plays a crucial role. More recently, it has been sug- 
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gested that the formation of partially ordered anisotropic 
domains and microphase structures may be resp onsible 
for physical gelation in these systems [HI, [H, 13 ■ 

In a recent letter [l3fl. we have introduced a <fi 4 model 
with competition between a short-range attraction, de- 
scribed by the Ginzburg-Landau Hamiltonian, and a 
long-range screened repulsion, described by a Yukawa 
potential |15j. Albeit schematically, the model contains 
the essential features of the effective interaction potential 
among charged colloids in polymeric solutions and sheds 
new light on the structural properties of these systems. 
Here we present a detailed and extensive analytical study 
of this model, and a full derivation and explanation of the 
results previously reported briefly in [131 ]. 

We show that depending on the control parameters, 
there is a region of the phase diagram where usual 
phase separation between a colloid-rich and a colloid- 
poor phase takes place. Conversely, as the screening 
length and/or the strength of the repulsion exceeds a 
threshold value, phase separation is prevented. In this 
case, at moderately high temperature the competition 
between attraction and repulsion has the effect to pro- 
duce microphase structures, i.e., partially ordered mod- 
ulated domains which in the terminology of particle sys- 
tems, correspond to the cluster phase. These modulated 
structures are the precursors of a first order transition 
to an equilibrium lamellar phase found at lower temper- 
atures. By using a replica approach for systems without 
quenched disorder [n| an d by employing the Self- 
Consistent Screening Approximation (SCSA) [H|,[l!|, we 
also show the presence of a glass transition line in the low 
temperature region, once the first order transition to the 
lamellar phase is avoided. The mechanism responsible 
for the glass transition in this case turns out to be com- 
pletely different from that of molecular glass-formers: as 
a matter of fact, we find that glassiness is not due to 
the presence of a hard-core interaction, which is totally 
absent in our model; it is instead due to the formation 
of the microphase structures which order up to the size 
of correlation length. The geometric frustration, arising 
from arranging such modulated structures in a disordered 
fashion, leads to a complex free energy landscape and, 
consequently, to a dynamical slowing down. 

Our results suggest that the cluster phase observed in 
colloidal suspensions should be followed, upon decreasing 
the temperature (or increasing the volume fraction), by 
an equilibrium periodic phase (a columnar or a lamellar 
phase, depending on the volume fraction). If, instead, 
this ordered phase is avoided, a structural arrest, corre- 
sponding to the gel phase observed in the experiments 
and in numerical simulations, should eventually occur. 
The existence of ordered phases in colloidal suspensions 
and the fact that the transition to the gel phase could 
occur in a metastable liquid, are novel predictions, which 
have never been considered before. Recently, the pres- 
ence of such ordered phases in atomistic model systems 
of charged colloids interacting via the DLVO potential, 
has been unambiguously shown both in three [ill . [lij and 



in two [111. Il2j dimensions. 

The paper is organized as follows: in Sec. [H] we de- 
scribe the model; in Sec. IIIII we study the model within 
the self-consistent Hartree approximation; in Sec. IIVI we 
analyze the dynamics of the system in the spherical limit; 
in Sec.[V]we study the glass transition; in Sec. lVII we com- 
pute the pair correlation function of the system in the 
real space, showing the emergence of competing length 
scales and emphasize the prominent role of microphase 
structures; in Sec. IVHl we summarize the results found: 
we discuss the resulting phase diagram and stress the 
connections with charged colloids. Some details of the 
calculations are reported in Apps. [AUDI 



II. THE MODEL 

We consider the standard three dimensional 4 field- 
theory with the addition of a repulsive long-range poten- 
tial nl: 



H[tj>] = 



d 3 . 



/(</>) 



\ (V^>(x)) 2 



+ H LR [<I>], (1) 



where 0(x) is the scalar order parameter field, related 
to the concentration of colloidal particles. The local free 
energy has the usual Ginzburg-Landau form: 



m = Y (0(x)) 2 + | (0( X )) 4 . 



(2) 



The coefficient ro is a temperature dependent mass pro- 
portional to the deviation ro oc T — T^ IF from the mean 
field transition temperature in absence of frustration. 
The long-range repulsive interaction is described by a 
Yukawa potential: 



W 



d 3 xdV 6 ' X _ X x % (x)^(x'). (3) 



The local term favours the formation of a uniform con- 
densed phase, while the long-range term energetically 
frustrates this condensation. In the momentum space 
the Hamiltonian of Eq. ([I]) reads: 



d 3 k 



r + k 2 + 



inW 



A- 2 + 



k 2 



0k (, 



(4) 



gV f d 3 ki d 3 k 2 d 3 k 3 



(2tt) 3 (2tt) 3 (2tt) 3 



The model has been also studied in [20( in the context 
of microemulsion. The parameters W and A are, respec- 
tively, the strength and the range of the repulsive po- 
tential. For W = we obtain the canonical short-range 
ferromagnet. Interestingly, for A — ► oo we recover the 
case of Coulomb interaction. The latter model was first 
introduced in Ref. [2l| . in the context of cuprate sys- 
tems. Then, it has been further studied in many other 
papers (see, for instance, Refs. [H, [H, [H, H, [H [13] ) 
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FIG. 1: (color online) Main frame: Phase diagram of the 
model obtained within the self-consistent Hartree approxima- 
tion, in the frustration (47rW / )-temperature (T) plane, for a 
fixed range of the repulsive potential, A = 2 (and g = 1, 
ro = —1). The temperature is rescaled with respect to the 
critical temperature of the unfrustrated model, T C (W = 0) = 
27r 2 /3A. The (red) continuous curve, T C (W, A), corresponds 
to the second order phase transition from the paramagnetic 
to the ferromagnetic phase (with the usual Hartree critical ex- 
ponents, v = 1 and 7 = 2), whereas the (blue) dashed curve 
corresponds to the fluctuation-induced first order transition 
from the paramagnetic to the lamellar phase, Tl(W, A). W c 
is the critical threshold of the repulsion strength, Eq. 
5are set in such a way that the elementary length of a lat- 
tice spacing, is Inset: System phase diagram in the range 
(A)-frustration (W) plane for a fixed temperature T = 0.1, 
showing the relative position of the different phases. 



where it has been used to describe the phenomenology 
of a wide variety of systems, where competing interac- 
tions on different length scales stabilize pattern forma- 
tion and the creation of spatial inhomogeneities (for a 
review see [28]). These systems include magnetic systems 
and dipolar fluids characterized by long-range Coulom- 
bic interactions , mixtures of block copolymers [3(| , 
water-oil-surfactant mixtures [3l| and do ped Mott insu- 
lator, including high T c superconductors [321 ]. As a con- 
sequence, our model allows to describe and to interpret 
in an unified fashion the phenomenology of a broad range 
of different physical systems. 



consistently [33j. For a one component order parameter 
there are six ways of choosing the two factor of </> to be 
paired in (cf> 2 ) among the four factor of </> in the quar- 
tic term <?</> 4 /4. Hence, the self-consistent Hartree ap- 
proximation consists in substituting the term gcj. r/4 with 
3<?(0 2 )</> 2 /2. With this substitution, the Hamiltonian be- 
comes quadratic: 



Z = Vfa e 



. f 3v 



( 5 ) 

It is then possible to compute Z and evaluate the corre- 
lation function: 

G(k) = (0 k 0-k)-(0k)(0-k>. (6) 
In the paramagnetic phase, ((/>k) = 0, one obtains: 

T 



G(k) 



ro 



T 



7,2 I 4:7rW 



A-kW 



(7) 



where the renormalized mass term, r, has been defined 
as: 



r = r + ig{4> 2 ). 



We now recall that: 



G(x,x) = f 

J\k\<A 



d^k 

(2tt) 3 



G(k), 



(8) 



(9) 



where A is the ultraviolet cutoff. As a result, Eqs. © 
and © yield: 



r = r + 3.9 



|k|<A 



d^k 

(2tt) 3 



k 2 



A- 2 + fc 2 



(10) 



which provides a self-consistent equation for the renor- 
malized mass r. It can be shown that, apart from the 
factor 3, this approximation becomes exact for an N- 
component order parameter in the limit N — > oo 33] . 

In order to solve the self-consistent equation, Eq. (|10[) . 
it is convenient to define: 



1 



4ttA 4 ' 



(11) 



and distinguish between two cases: W < W c and W > 



III. PARAMAGNETIC, FERROMAGNETIC 
AND LAMELLAR PHASES WITHIN THE 
HARTREE APPROXIMATION 

In this section we solve the model, Eq. |T]), within the 
self-consistent Hartree approximation. This approxima- 
tion amounts in replacing one factor tfi 2 in the quartic 
term by its average, namely (4> 2 ), to be determined self- 



A. W < W c 

For W < W c the correlation function, G(k), behaves as 
in the standard unfrustrated case: it is a monotonically 
decreasing function of |k|, with a maximum in k = 0, 
where it equals: 

G(k = 0) = -. (12) 
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Therefore, as r — > —4irWA 2 , the susceptibility, x = 
G(k = 0)/T, diverges. This corresponds to a usual 
second order phase transition towards a ferromagnetic 
phase. From Eq. (|10[) . one can compute the critical tem- 
perature, T C (W, A), given by the following relation: 



3gT c 

2tt 2 



dk 



1 + fc 2 



1 



ill 



= -r - 4ttPFA 2 . (13) 



It is possible to show that the usual Hartree critical ex- 
ponents are found, i.e., v = 1 and 7 = 2 (see App. [A"]) . 
In the terminology of colloidal suspensions, this means 
that for T <T C the system undergoes a phase separation 
between a colloid-rich and a colloid-poor phase. The de- 
pendence of the critical temperature, T C (W, A), upon W 
(for A = 2 and ro = —1) is plotted in the main frame of 
Fig. [1] (red continuous line in Fig. [T]), showing that the 
only effect of the repulsive interaction for W < W c is to 
decrease the numerical value of T c : for W = we recover 
the usual critical temperature of the unfrustrated sys- 
tem in the Hartree approximation, T C (W = 0) = — 2 ^ g ^° 
(notice that tq < 0), whereas T c vanishes as W — > W c . 



B. W > W c 

Conversely, if W > W c , the correlation function has a 
maximum around a finite value of the wave vector, k m '■ 



1/2 



VlirW - A~ 



(14) 

As we shall discuss later on, the peak in the correlator 
indicates that the system establishes microphase struc- 
tures with inverse domain size given by k m , which char- 
acterize the incipient periodic order. Such modulated 
structures are the analog of the cluster phase observed 
in colloidal systems characterized, also in this case, by a 
peak in the structure factor around a characteristic wave 
vector correspondingto the inverse of the typical size of 
the clusters @, 0, B 1, (3 ■ 

In k = k m , the correlation function is given by 



G(|k| = k n 



+ 2V^W - A- 2 



(15) 



Therefore, G(|k| = k m ) and, consequently, the fluctua- 
tions, x(|k| = k m ), diverge for r — > r sp = A~ 2 — 2\J A_-nW . 
As a result, also the integral of the right hand side of 
Eq. (|10p diverges, and the self-consistent Hartree rela- 
tion can be satisfied only at T = 0. This implies that 
for W > W c the paramagnetic phase is stable at all fi- 
nite temperatures, its spinodal line, where the suscepti- 
bility for |k| = k m diverges, x(fc m ) — * 00, is located at 
T = 0. Therefore above W c there is no phase separation. 
This result is quite important for designing new materi- 
als as well as in the experimental and numerical study 
of colloidal systems, where it is crucial to distinguish the 



slowing down due to colloidal gelation from that due to 
kinetic of phase separation. 

Nevertheless, the system still undergoes a first order 
transition to a lamellar phase when the temperature is 
lowered below Tl(W, A) (blue dashed line in Fig.Q}. Such 
first order phase transition is induced by the fluctuations, 
as first discussed by Brazovskii for a related model 1341 
and as also found in the Coluombic case (A —> 00) [271 ] . 
The transition temperature Tl(W, A) can be determined 
on equating the free energy of the paramagnetic phase to 
that of the lamellar phase [27|, as explained in App. [B] 

The lamellar phase is characterized by a periodic vari- 
ation of the order parameter with wave-length l m = 
2-Kkz}: 



t>k) = m (5 (k - k m ) + 6 (k + k m )) , 



(16) 



where k m is given in Eq. (|14p . Note that as W approaches 
W c from above, according to Eq. (Tj"4")) . k m vanishes. This 
signals the fact that the size of the stripes diverges on 
the boundary between the lamellar and the ferromagnetic 
phase. 

Within the fluctuation approach described in App. [Bl 
it is also possible to study the stability of other kinds of 
ordered phases, such as columnar or periodically ordered 
cluster phases, which occur in numerical simulations [T3 , 
Hil ]. In fact, it turns out that these phases can be stable 
besides the lamellar phase, at lower volume fractions, i.e., 
for (0k=o) < 0. However, in this paper we only focus on 
the case (</>k=o) = 0, (i.e., J d 3 x(0(x)) = 0) where the 
only periodically ordered stable phase is the lamellar one. 



C. Phase diagram within the Hartree 
approximation 

The results found within the self-consistent Hartree ap- 
proximation are summarized in Fig. [TJ where the system 
phase diagram in the frustration (47rVF)-temperature (T) 
plane at fixed screening length A = 2 (main frame) and 
in the screening length (A)-frustration (W) plane at fixed 
temperature T = 0.1 (inset), is plotted, showing the rela- 
tive positions of the paramagnetic (P), the ferromagnetic 
(F) and the lamellar (L) phases. 

According to Eq. (fTl"]) . by changing A, W c changes and, 
hence, the phase diagram modifies. For instance, in the 
Coulomb limit, A — * 00, we have that W c — > 0. As a 
consequence, the ferromagnetic phase reduces to the axis 
W = 0, for T < T C (W = 0). Conversely, for A -> we 
have that W c — > 00 and the lamellar phase disappears. 
It is interesting to observe that the dependence of the 
threshold W c on the screening length A is exactly the 
same of that found in numerical simulations of a model 
system for charged colloids interacting via a DLVO po- 
tential 0: the authors find that for W < W c the sys- 
tem tends to phase separate, whereas for W > W c the 
preferred structures are elongated (and eventually quasi 
one-dimensional) finite clusters. Physically, such thresh- 
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old value, W c , can be interpreted expanding the Hamil- 
tonian, Eq. (0]), up to the lowest order in k 2 : 



„, rn V f d 3 k 



AttWX 



(17) 



gV f d J kx d J k 2 d A k 3 



(2tt) 3 (2tt) 3 (2tt) 3 ^l^VKaV-Ki-Ma-K.- 

At W = Wc the coefficient of the k 2 term vanishes, im- 
plying that there is no energetic cost associated to the 
creation of interfaces between high density and low den- 
sity regions, leading to pattern formation and to the cre- 
ation of spatial inhomogeneities. 



IV. DYNAMICS IN THE SPHERICAL LIMIT 

In this section we study the dynamics of the model in 
the case of an iV-component order parameter in the limit 
N — > oo, after an instantaneous quench from high tem- 
perature. The spherical approximation is equivalent to 
the Hartree approximation described in the previous sec- 
tion (apart from the factor 3 in Eq. ([TO]) ), and allows to 
study analytically the dynamics. We solve the time de- 
pendent Cahan-Hilliard equation [35| . following the ap- 
proach taken in Refs. [36| . We obtain an exact analytical 
solution for the evolution of the time dependent struc- 
ture factor of the system, G(k,t), and for the peak po- 
sition, k max (t). In agreement with the results presented 
in the previous section, we find that as W < W c , the 
usual spinodal decomposition occurs, and the maximum 
of the structure factor approaches zero as t- 1 / 4 for large 
times. This corresponds to phase separation between a 
high density colloid-rich phase and low density colloid- 
poor one. On the other hand, for W > W c , provided 
that we are on the spinodal line of the lamellar phase 
(T = 0), the peak position approaches a non-zero value, 
k m , given in Eq. (fTi]) . and the structure factor asymptot- 
ically approach a 5- function peaked around k m , signaling 
a lamellar order. However, since the paramagnetic phase 
is stable at all finite temepratures, for T > we find that 
the lamellar phase starts to form locally but then fades 
away after long time and the system stays homogeneous: 
kmax(t) approaches a limiting value where the structure 
factor vanishes exponentially. 

Let us consider the Cahan-Hilliard equation [35] : 



d0(x,<) 

at 



= MV 



; ^[0(x,t)] 

^(x,t) ■ 



(18) 



where M is related to the mobility, and TL is given in 
Eq. Ilj . For an A-component order parameter the above 
equation reads: 



dt 



= MV 1 



N 



r + If 

V 0=1 



V 2 U Q (x,f) 



w 



, e -|*-x'|A0 a ( x ' )t ) 



(19) 



In the limit N — > co we recover the spherical model 

replacing ^E^=i^( x .*) b y (^a( x >*))- Here rep- 
resents an ensemble average over the initial configura- 
tions. By assuming the translational invariance of the 
pair correlation function G(x, x',i) = (</>(x, t)0(x', t)) = 
G(|x-x'|,t), the quantity (0 2 (x,i)) = G(0,i) = S{t) 
becomes independent on the position x. Thus, dropping 
the index a, taking the Fourier transform over the space 
of Eq. (fT9"|) and multiplying both sides by 4> a {— k, t) we 
obtain the equation of motion for the time dependent 
structure factor, G(k, t) = (0 a (k, t)<j) a {— k, t)): 



ggCM) 
dt 



-Mk 2 r + gS(t) + k 2 + — G(k,t). 

A T - h 

(20) 

By replacing 2/3=1 4^ ( x ) t) by S(t) we have, in fact, 
enabled a linearization of Eq. (fTO]) . "preaveraging" the 
nonlinear term. Now Eq. ([20]) can be readily integrated 
to obtain: 



G(k,i) = G(k,0)e 



-Mk 



[Q(t)+k- 



t+- 



where Q(t) is defined by 

Q(t) = f dt' [gS(t>) + r] 
Jo 



(21) 



(22) 



To complete the solution, S(t) must be computed self- 
consistently. As shown in App. [Cl one can obtain 
an expression for the time dependent structure factor, 
Eq. (|CJ7|) . and an equation for the evolution of the peak 
position: 



exp 



Alt 1 



= P x 



(A- 2 
1 - 



k 2 

max 



4ttW\~ 2 



1/2 



(23) 



(A- 



h 2 

max 



(*)) 3 



where the function P is given in Eq. (|C5jl . Taking the 
logarithm of both sides of Eq. (j2"3")) and dividing by Mt we 
obtain an asymptotic expression for k max (t), Eq. l|C8p . 
which is readily analyzed in the limit t — > co. Careful 
analysis show that depending on the values of W and A, 
different solutions of this equation are found. In partic- 
ular, one needs again to distinguish between the cases 
W < W r and W > W r . 



by: 



A. W < Wc 

For W <W C the solution of Eqs. (|23]) and ([C8]) is given 
fcoo = lim k max {t) = 0, (24) 
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corresponding to the usual spinodal decomposition. Us- 
ing Eq. (|C8p we find that that the usual N — > oo domain 
growth exponent of the canonical short-range ferromag- 
net with conserved order parameter is found: 

k m ax{t) ~ t / . (25) 

Obviously, this solution can be obtained provided that we 
are below the critical point, r < -4ttW\1 [T < T C (W)}, 
otherwise no spinodal decomposition takes place and the 
system remains in the paramagnetic phase. 

B. W > W c 

For W >W C the solution of Eqs. and (|C8)l is given 
by: 

kl = - A" 2 = (26) 

where the wave vector fcoo coincides with that defined in 
Eq. (fT4|) . In this case, the peak position reaches a nonzero 
steady state value for large times and the structure factor 
asymptotically approaches a ^-function peaked around a 
finite value of the wave vector k m . This situation cor- 
responds to the lamellar phase and k m is the inverse 
domain size. In particular, for earlier times (such that 
AttW < k max (t)), using Eq. JC8|, we have that: 

fcmax(i) ~ UlntJ t" 1/4 . (27) 

Thus, apart from the logarithmic factor, in this regime 
k m ax(t) behaves as the standard solution of the unfrus- 
trated short-range ferromagnet in the large N limit in 
the case of usual spinodal decomposition. However, at 
later times, k ma x(t) saturates at k m as a power law: 

fcL*(*)~*4 + C*- 1/2 . (28) 

According to Eq. (|C7[) , the structure factor G(k, t) can 
be approximated by a Gaussian centred about k rn and 
with width 

Ak~- [ -^—^ T , (29) 

which in the limit t — > oo approaches a 5-function centred 
in k m , as stated before. 

However, this situation can occur only provided that 
we are on the spinodal line of the lamellar phase, located 
at T = 0. As discussed in the previous section, the ex- 
pression of the spinodal line reads (both in the Hartree 
approximation and in the large N limit): 

r sp (W,X) = -2V4vrVK + A~ 2 . (30) 

Indeed, as r > r sp (W, A) (i.e., T > 0), from Eq. (|C8|) we 
have that k max (t) approaches exponentially a limiting 
value, £4 = -{r + A~ 2 )/2: 

e -M(W-W sp )t 
k m ax(t) = ^oo + C x 7yt= ' (31) 




FIG. 2: (color online) Main frame: Complexity, E, as a 
function of T/Tk for 4^14^ = 0.5 (red dashed curve) and 
4ttW = 0.08 (black continuous curve), for a fixed value of the 
range of the repulsive potential, A = 2, and for g = 1 and 
ro = — 1. At high temperature (T > Td) the system is in the 
paramagnetic phase and the complexity is zero. At T — Td 
the complexity discontinuously jumps to a finite value, signal- 
ing the emergence of a complex free energy landscape. The 
complexity decreases as the temperature is decreased and van- 
ishes at Tk where the thermodynamic transition to a 1RSB 
glassy phase takes place. Inset: Complexity at the dynami- 
cal transition temperature, T,(Td), plotted as a function of the 
strength of the repulsive potential, 4nW. As W decreases, 
£(Td) decreases. At W = W c , the glassy phase disappears. 

where W sp = (A~ 2 -rf/lQir. Consequently, P, Eq. (|C5|) . 
asymptotically vanishes in an exponential fashion, signal- 
ing that the periodic pattern fades. At the spinodal line 
we have that W — W sp , which gives back Eq. I|28p. 

In conclusion, the analysis of the dynamics in the 
spherical limit indicates that for W < W c the system un- 
dergoes a phase separation between a colloid-rich phase 
and a colloid-poor one. On the other hand, for W > W c 
and T > the system starts to form microphase struc- 
tures and periodically ordered modulated domains, as 
signaled by a peak in the structure factor which grows 
in time and whose position seems to approach k m . How- 
ever, since the paramagnetic phase is stable at all finite 
temperatures, the microphase structures have a finite 
lifetime. The lamellar order is fully established only at 
T = 0. 



V. THE GLASS TRANSITION 

Since the paramagnetic phase is stable at all finite tem- 
peratures, the first order transition to the lamellar phase 
can be kinetically avoided. It is then possible to super- 
cool the system in the (meta-stable) homogeneous phase 
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below the first order transition line, Tl(W,X), as much 
as one can obtain supercooled liquids by performing fast 
enough coolings below the melting temperature avoiding 
crystallization. Experimental [37j and numerical [38l |39| 
results show that, indeed, these systems exhibit long-time 
relaxations similar to that observed in glasses. Recent re- 
sults have clearly shown the presence of a glas s transition 
in the Coulomb case (A — > oo) [25|, [26j, 1271. 140]. 



A. General formalism 

In order to study the glass transition in our model, we 
use a replica approach introduced in Ref. [l6| and then 
developed in Ref. [l7|, formulated to deal with system 
without quenched disorder, which allows to compute the 
complexity, E. 

Our aim is to derive the free energy landscape of the 
model. The equilibrium free energy, defined as F — 
— TlniT, is relevant only if the system is able to explore 
the entire phase space. This is not the case in the glassy 
phase, where the system is frozen in metastable states. 
In order to scan the locally stable field configurations, we 
introduce an appropriate symmetry breaking field ?/>(x), 
and compute the following partition function [161 ] : 



V$ exp -iiUW\ - 7T / d J x [V>(x) - 0(x)] 



where u denotes the strength of the coupling. Since u > 
0, the introduction of the pinning field ip(x) forces the 
order parameter </>(x) to assume configurations close to 
-0(x). Therefore, the free energy 



M =-Tln%] 



(33) 



will be low if tp{x) equals to configurations which locally 
minimize H[<fi]. Thus, in order to scan all metastable 
states, we have to sample all configurations of the field 
V>, weighted with exp(— f3f[ip]): 



F = lim 

u->0+ 



fThfifty] exp (-/?/>]) 
JP^exp (-/?/>] 



(34) 



Therefore, F is a weighted average of the free energy 
in the various metastable states. If there are only a 
few local minima (i.e., a non-extensive number of local 
minima) of free energy / equal to F, we have that F 
equals the true free energy F of the system. However, in 
case of the emergence of an exponentially large number 
of metastable states with large barriers between them, 
a nontrivial contribution arises from the above integral 
even in the limit u — > + and F differs from F. This 
allows to identify the configurational entropy E [H, [13]: 



In order to get an explicit expression for E we introduce 
replicas. The replicated free energy reads: 



F(m) = - lim —In [ Vi\) [z\£\\ 
«^o+ m J V / 

from which, evidently, F can be obtained as 

dmF(m) 



dm 



F 

and hence, according to Eq. ([55] 
1 dF(m) 



E = 



T dm 



(36) 



(37) 



(38) 



Using Eq. (f52j) and integrating over ifi, we get: 

„ rn / m 

Z(m) = lim + / Y[V(j) a exp - (3^H[(j) a ] 

U a—1 \ a=l 

J2 |d 3 x^(x)/(x)j. (39) 



u 

2m 



The partition function of Eq. (|39|) is formally equivalent 
to that of a system in a quenched random field analyzed 
by means of the replica trick, such as the random field 
Ising model (RFIM). The only difference is that in this 
case the limit m — > 1 must be taken. We can thus use the 
techniques developed to deal with such systems [l9( . The 
matrix correlation function of the problem of Eq. 
obeys the following Dyson equation: 



G- 6 x (k) = Go x (k)5 o6 



E Qb (k) , 

m 



(40) 



F = F - TE. 



(35) 



where G ~ 1 (k) is the propagator of the free unreplicated 
theory, Eq. (J7J), and E a t(k) is the self-energy in the 
replica space. If we find that, due to the ergodicity- 
breaking coupling constant u, E Q b(k) has finite off- 
diagonals elements, we can conclude that there must be 
an energy landscape sensitive to that infinitesimal per- 
turbation, leading to a glassy dynamics. On the other 
hand, if the self-energy is diagonal, the dynamics is er- 
godic and the system is in the liquid phase. 



B. The Self-Consistent Screening Approximation 

In the following we use the self-consistent screening 
approximation (SCSA) [H, [l9|, [2^, [2(| , which consists in 
introducing a A^-component version of the model with 
4> = (0i, . . . , 4>n) and coupling constant g/N and sum- 
ming self-consistently all the diagrams of order 1/N (see 
Fig- 13 ■ This approximation is exact up to order 1/N. At 
the end of calculations we will consider the scalar case, 
N = 1 41]. Since the attractive coupling between repli- 
cas is symmetric with respect to the replica index, we 
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FIG. 3: Main frame: Momentum dependence of the corre- 
lation function, G(k), for 4ttW = 0.2 and A = 2 at T = T d , 
showing that it is peaked around the typical modulation wave 
vector km with width given by the inverse of the correla- 
tion length, signaling the fact that microphase structures es- 
tablish over a finite range f. Inset: Momentum dependence 
of the non ergodicity parameter, /k, for the same values of 
W, X and T. 

assume the following structure of the matrix correlation 
function: 

G afc (k) = [0(k) - JF(k)] S ab + T(k), (41) 

i.e., with diagonal elements <?(k) and off-diagonal ele- 
ments .F(k). It can be shown that for system with 
quenched disorder (such as the RFIM), this ansatz turns 
out to be equivalent to the one-step replica symmetry 
breaking ansatz (1RSB). While the diagonal correlator 
can be interpreted as the usual (one time) correlation 
function, i.e., 

TG(k) = (fat-n), (42) 

the off-diagonal term can be interpreted as measuring the 
long-time correlations: 

TJF(k)= hm(0 k (^_ k (O)}. (43) 

t — >OG 

Hence, .F(k) vanishes in the paramagnetic phase while is 
finite in the glassy one. 

As shown in Ref. [26[ and in App. [D] within the 
SCSA the expressions of the correlators can be found by 
solving numerically a set of coupled integral equations, 
Eqs. |D8)) - (|DT6l) . After that, according to Eqs. jDl~9l) 
and (1D20|) we are able to compute the complexity, E. 

C. Results 



ity for different values of the strength of the repulsion, 
W, and of the temperature T. The model undergoes 
a glass transition of the same nature of that found in 
discontinuous spin glasses 

[HEjllii. In Fi &® the be- 
haviour of the configurational entropy is plotted for two 

different values of W as a function of the temperature. 
For a given value of the strength of the repulsion W , at 
high temperature, ^"(k) vanishes, leading to a vanishing 
complexity, corresponding to the fact that the system 
is in the liquid phase. At T = T d {W) the configura- 
tional entropy jumps discontinuously to a finite value, 
signaling the emergence of an exponentially large num- 
ber of metastable states. At this temperature a glassy 
dynamics sets in, corresponding to a non-zero value of 
the long time correlation function .F(k). The complexity 
decreases as the temperature is decreased and vanishes 
at Tk{W). At this temperature, called the Kauzmann 
temperature, an ideal thermodynamic glass transition to 
a 1RSB glassy phase takes place. Note that, in agree- 
ment with the results presented in the previous section, 
the glass transition disappears for N — > oo, when the cor- 
rections due to the diagrams of order 1 /N are not taken 
into account. 

In order to characterize the nature of the glass transi- 
tion, we examine the properties of the correlation func- 
tions In the main frame of Fig. [3] the diagonal part 
of the correlator, <?(k), is plotted as a function of the 
wave vector at the dynamical transition temperature Td, 
showing that the correlation function is clearly peaked 
around a maximum located in k m , whose value is given 
in Eq. (fT4|) . As shown in the next section, the width 
of such maximum, allows to identify the inverse of 
the correlation length. The shape of <?(k) indicates that, 
although no periodic order occurs ((</>k m ) = 0), a lamel- 
lar structure of wave length l m = 2irk^ n 1 over a finite 
range £ is formed. These microphase structures are the 
analog of the cluster phase observed in colloidal systems. 
The system forms a mosaic of such modulated structures, 
periodically ordered over length scales smaller than the 
correlation length £, and randomly assembled in a dis- 
ordered fashion over larger length scales. The defects 
and the imperfections in the perfect stripe arrangements 
give rise to the tails in the correlation function. The 
lower is the temperature, the more pronounced is the 
peak of £7(k) around the maximum in k m , signaling the 
fact that £ increases as the temperature is decreased. At 
T = Td, where £ > 2l m , these microphase structures es- 
tablish over a length larger than their modulation length 
and becomes frozen. The glass transition arises from the 
fact that there are many possible configurations to ar- 
range such modulated structures in a disordered fashion, 
leading to the emergence of metastable states. 

The characteristic wave length, fc~ , dominates also 
the dynamics as indicated by the momentum dependence 
of the non ergodicity parameter 



We have fixed the range of the repulsive potential to 
A = 2 and we have studied the behaviour of the complex- 
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FIG. 4: (color online) Main frame: Dynamical (blue curve 
and triangles) and thermodynamical (red curve and dia- 
monds) transition temperature to the glassy phase, as a func- 
tion of the strength of the repulsive potential, 4nW, for A = 2. 
Inset: Values of the renormalized mass, rd and rjf, corre- 
sponding, respectively, to the dynamical and the thermody- 
namical glass transition temperatures, as a function of 4nW. 
In the figure the difference between rd(W) (red curve and tri- 
angles) and tk{W) (blue curve and diamonds) and the mini- 
mum value of the renormalized mass, r ap (W), is plotted. The 
figure shows that as W — > W c , both rd and tk approach 



,(W), signaling that T d , K {W ) 



0. 



plotted in the inset of Fig. [3] at T = Td- Basically, /k 
represents the height of the plateau reached by the cor- 
relation function of the density fluctuations of wave vec- 
tor k after infinite time. From Eq. (|44[) one immediately 
observes that as = (T > Td), also the non ergod- 

icity parameter, /k, vanishes: the system is in an ergodic 
phase, since all the correlation functions relax to zero. 
On the other hand, at T = Td the ergodicity is broken. 
The presence of a maximum of /k located in |k| = k m 
signals the fact that structural arrest is more pronounced 
over length scales of order l m . The width of the maxi- 
mum of the non ergodicity parameter, T _1 , can be tied 
to the emergence of a new length scales. Indeed, the 
height of the plateau decays as ||k| — k m \ > T _1 , and T 
represents the typical length scale over which defects and 
imperfections in the periodically modulated pattern are 
allowed to wander and to diffuse over long times [25|, H|| . 

In Fig. 0] the dynamical transition temperature Td(W) 
and the ideal one Tk{W) are plotted as a function of 
the frustration 4irW, for a fixed value of the range of 
the repulsive potential, A = 2. It is also interesting to 
study the behaviour of rd(W) and r# (W), i.e., the values 
of the renormalized mass r corresponding, respectively, 
to the dynamical and to the thermodynamical transi- 
tion temperature. More precisely, in the inset of Fig. 31 
rd(W) — r sp (W) and r K (W) — r sp (W) are plotted, where 



FIG. 5: (color online) Temperature dependence of the char- 
acteristic length scales. At high temperatures, T > T* , the 
system behaves as a fluid of charges of linear size £ (red curve) 
and screening length Id (black curve). At T* this description 
breaks down and the competition between attraction and re- 
pulsion produces microphase structures with interstripe dis- 
tance l m (blue curve), established over a finite range £. The 
modulation length, l m , decreases as the temperature is de- 
creased and at low temperature equals 27rfc~ 1 . Conversely, 
the correlation length, £, increases as the temperature is de- 
creased and diverges at T — » at the spinodal line of the 
homogeneous phase. Classiness emerges when £ > 2l m . On 
the contrary, the modulated structures progressively fade as 



fsp{W) is the minimum value of the renormalized mass 
(corresponding to T = 0), whose expression is given in 
Eq. ((30)) . The figure clearly shows that the glassy phase 
disappears as W — > W c , where we find that rd t K{W c ) — > 
r sp (W c ), i.e., T d {W c ) = T K (W C ) = 0. (Nevertheless, we 
have that Iim lf _ >w + T d , K {W) ^ T d , K (W c ) = 0). These 
results can be interpreted if one observes that, according 
to Eq. fl4"]) . the size of the modulated structures increases 
as W is decreased and diverges as W — > W c . In fact, for 
purely geometrical reasons, the larger is the modulation 
length, the smaller is the number of metastable configura- 
tions to arrange the microphase structures in a disordered 
fashion. As a consequence, the number of metastable 
states (i.e., the complexity) vanishes as W — * W c and 
the glassy phase disappears. This is clearly outlined by 
the inset of Fig. ^ which shows that the complexity com- 
puted at the dynamical transition temperature tends to 
zero as W — > W c . 
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VI. MODULATED STRUCTURES AND 
CHARACTERISTIC LENGTH SCALES 

In this section we analyze more accurately the nature 
of the microphase structures which are formed due to 
the competition between attraction and repulsion on dif- 
ferent length scales. This investigation provides many 
insights on the nature of the cluster phase observed in 
colloidal suspensions and on the physical mechanisms in- 
ducing the structural arrest in these systems. 

Let us consider the case W > W c , where the correla- 
tion function, G(k), is peaked around the typical mod- 
ulation wave vector, k m . In the following we compute 
the correlation function in real space (within the self- 
consistent Hartree approximation, or, equivalently, in the 
spherical limit) and we show that different length scales 
emerge in the system [2^, [42| ■ The pair correlator in real 
space is determined by computing the Fourier transform 
of G(k), given in Eq. 0: 



G(X) = (2^/ d3kG(k)elkX 



(45) 



T f „ k(\- 2 + k 2 ) sin(fcr) 
are ■ 



2tt 2 |x|7 (k 2 + ( 2 )(k 2 + k 2 ) 



where 



C 2 K 2 = (r + \- 2 )±V(r-\- 2 ) 2 ~l6nW ^ 

The above equation marks a crossover temperature T* , 
corresponding to the temperature at which the renormal- 
ized mass equals: 



r* = A" 2 + 2V4ttW. 



(47) 



A. High temperature: 



r > r 



For high enough temperatures, r > r* (i.e., T > T*), 
the argument of the square root of the right hand side of 
Eq. (|4*o| is positive. Consequently ( 2 and k 2 assume real 
values and the integral, Eq. (|45|) . can be easily evaluated 
by applying the residue theorem to the poles lying on the 
imaginary axis at k = ±i£, ±i«, leading to 

r , :] _ T[(( 2 -X- 2 )e-^ + (X- 2 - K 2 )e-^] 

G(X) " in\K\(C 2 ~K 2 ) ' (48) 

The above expression clearly indicates the emergence of 
two characteristic length scales, £ = |3?{C}| _1 and Id = 
|9?{ft}| . The former plays the role of the correlation 
length of the canonical short-range ferromagnet and at 
high temperature behaves as: 



-1/2 



(49) 



The other length Id, instead, is a renormalized screening 
length, and for high temperature tends to 



the range of the repulsive potential. Therefore, the sys- 
tem at high temperatures, T > T* , behaves as a fluid 
of charges of linear size £ and screening length Id ■ The 
temperature dependence of £ and Id is shown in Fig. [5j 



B. Low temperature: r < r* 

Conversely, for r < r* (i.e., T < T*) the argument 
of the square root of Eq. (I46[) becomes negative. Conse- 
quently, £ and k 2 assume complex values. The analytic 
continuation of Eq. (|48|) to the low temperature region 
reads: 



G(x) 



T 



{2k\K2 — A 2 ) cos(k2|x| 



167t|x|kiK2 

+ [k\ -k 2 - \- 2 ) sin(K 2 |x|)| , (51) 



where k — k\ 4- in2- This expression implies that, al- 
though no periodic order occurs ((</>k m ) — 0), a lamellar 
structure of wave length l m over a finite range £ is formed. 
The modulation length is given by: 



l m = 2tt|k 2 | 1 = 4tt 



2 VrA" 2 + A-kW — r - A~ 2 



-1/2 



(52) 

whereas, according to Eq. (fSTj) . the correlation length 
reads: 



(r + A" 2 ) + 2^/r\- 2 + AirW 



-1/2 



h 



A, 



(50) 



and gives the range over which the microphase structures 
are formed. The temperature dependence of the charac- 
teristic length scales present in the system is shown in 
Fig. El The correlation length increases as the tempera- 
ture is decreased and diverges as T goes to zero, where 
the spinodal line of the paramagnetic phase is located. 
At low temperatures we have l m = 2irk^-, where k m is 
the characteristic wave vector defined in Eq. (I14|) . Con- 
versely, as T — > T*~, the modulation length l m diverges 
as (r* — r)~ x l 2 . Interestingly enough, it turns out that 
the glass transition emerges at temperatures at which 
£, > 2i m , as conjectured in Ref. [261 ] : the microphase 
structures form over a length larger than their modu- 
lation length and become frozen. The glass transition 
arises from the fact that there are many possible config- 
urations to arrange such modulated structures in a dis- 
ordered fashion, leading to a great number of metastable 
states. As the temperature is increased, the cluster phase 
continuously fades approximately at the temperature at 
which £ < l m : the microphase structures are ordered 
over a length scale smaller than their own modulation 
length and the system becomes homogeneous. Such a 
crossover temperature at which approximately the mod- 
ulated structures fades and the system is homogeneous 
corresponds to the black dotted line of Fig. [Sj 



11 




modulated 
structures 



W 0.1 
c 47CW 



0.2 



FIG. 6: System phase diagram in the frustration (4nW)- 
temperature (T) plane, for a fixed value of the screening 
length of the repulsive potential, A = 2, and for g = 1 and 
ro = — 1. The temperature is scaled with respect to the crit- 
ical temperature of the unfrustrated system, T C (W = 0) = 
— 27r 2 ro/3(?A. For W < W c a second order phase transition 
with the usual Hartree critical exponents from a paramagnetic 
phase (P) to a ferromagnetic (F) one is found, by lowering the 
temperature below T C (W, A) (red continuous curve). In the 
terminology of colloidal system this corresponds to a phase 
separation between a colloid rich and a colloid poor phase. 
For W > W c the phase separation is prevented and, instead, 
a first order transition to a lamellar phase (L), characterized 
by a periodic variation of the density, is found at Tl(W,X) 
(blue dashed line). Such a first order transition can be avoided 
and, in this case, at lower temperature a transition to a glassy 
phase (G) is found, characterized by a dynamical, Td(W,X) 
(green dotted curve), and by a thermodynamical, Tk{W, A) 
(magenta dashed-dotted curve) transition. The black dot- 
ted line represent a crossover temperature, corresponding to 
the temperature at which the correlation and the modulation 
lengths equal, £ ~ Z m , and the system establishes microphase 
structures, which are the analog of the clusters observed in 
colloidal systems. 



VII. PHASE DIAGRAM AND CONCLUDING 
REMARKS 



The results found in the previous sections are summa- 
rized in Fig. [6l where the complete phase diagram of the 
system is presented, showing the relative position of the 
different phases. The model, albeit schematically, retains 
the essential physics of charged colloids in polymeric solu- 
tions, where the competition between attraction and re- 
pulsion on different length scales plays a prominent role. 
The results found provide many insights on the equilib- 
rium phase diagram of colloidal suspensions, as well as 
on the physical mechanisms responsible for the formation 
of modulated structures and for the structural arrest. 

We point out that in our model, for W > W c , the com- 



petition between attraction and repulsion has the effect 
to produce microphase structures, which are the analog 
of the cluster phase observed in colloidal system, charac- 
terized by a peak in the structure factor around a wave 
vector corresponding to the inverse of the typical size of 
the clusters. The width of the peak is related to the cor- 
relation length, which represents the range over which 
such modulates structures are locally ordered. In our 
model the modulated structures are the precursors of an 
equilibrium lamellar phase occurring at lower tempera- 
ture. Therefore, according to this analogy, our results 
suggest that the cluster phase observed in colloidal sus- 
pensions might be the signature of the presence in the 
system phase diagram of an equilibrium lamellar phase, 
which is very often kinetically avoided on the experimen- 
tal time scales. Such an underlying equilibrium ordered 
phase may strongly affect the dynamics of the system in 
the low temperature region, close to the colloidal gela- 
tion. The presence of a stripe/columnar/lamellar phase 
in the phase diagram and the emergence of microphase 
structures has been clearly shown in recent numerical 
simulations of model systems of particles interacting via 
the DLVO pot ential both in two [TJ G3 and in three 
dimensions [14j , where the presence of first order transi- 
tions from a cluster phase to periodically ordered phases 
has been clearly detected. 

The nature of the glass transition found in our model 
gives many insights on the nature of the mechanisms in- 
ducing the phenomenon of colloidal gelation. We suggest 
that the structural arrest in colloidal suspensions might 
be due to the formation of microphase structures, which 
are ordered on small length scales, and are assembled 
together in a disordered fashion on larger length scales, 
thereby inducing a complex free energy landscape and, 
consequently, a complex and slow dynamics. 

Notice that in our model there is no geometric frus- 
tration due to the hard-core excluded volume, which is 
totally absent in the Hamiltonian. The frustration lies 
only in the competition between attractive and repul- 
sive interactions on different length scales. Therefore, 
our results seem to indicate that the physical mechanism 
inducing the structural arrest in charged colloids is to- 
tally different from that responsible for the jamming and 
the glass transition in molecular liquids at high volume 
fraction. 

In real experiments the addition of salt increases the 
number of ions in the solution, thereby decreasing the 
strength of the repulsive shoulder of the potential. Our 
results suggest that there is a critical value of the ions 
concentration (corresponding to W c ) above which no 
modulated structures are formed and the system under- 
goes, instead, a phase separation between a colloid-rich 
and a colloid-poor region. 

All these prediction could be experimentally and nu- 
merically checked. 
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APPENDIX A: CRITICAL EXPONENTS OF THE 
FERROMAGNETIC TRANSITION IN THE 
HARTREE APPROXIMATION 

In this appendix we compute the critical exponents of 
the second order ferromagnetic phase transition found 
within the Self-Consistent Hartree approximation for 
W < W c . 

As shown in Sec. lIIIl for W < W c , the correlation func- 
tion, G(k), has a maximum in |k| = 0, where it equals 
T/(r + AttWX 2 ). Therefore, as r -> r c = -4ttW\ 2 , the 
susceptibility, \ — = 0)/T diverges, corresponding 

to a second order critical point. Since the divergence oc- 
curs around |k| = 0, in order to determine the critical 
behaviour, we can expand the Hamiltonian of Eq. jlj up 
to the lowest orders in k 2 , as done in Eq. (fT7|) . neglect- 
ing the other irrelevant terms. According to Eq. (0), the 
susceptibility reads: 



and. thus: 



X 1 (r) = r + AttWX 2 = r - r c = t. 



(Al) 



Replacing r and r c with their expressions obtained via 
the self-consistent relation, Eq. ([10]) . and expanding the 
correlators up to 0(k 2 ), we get: 



3g_ 

2n 2 



dkk 2 



r + AttWX 2 + (1 - Ai:WX 4 )k 2 



{l-4irW\ 4 )k 2 _ 
Now, wc define the function I(t) as 

1 



J(r) 



A 

dfc — 

T- 



ck 2 



(A2) 



(A3) 



Recalling that r = r — r c = r + 4nW\ 2 = 1 — W/W c and 
that c = 1 - AkW\ a (note that c = for W = W c ), in 
terms of the function J(r) Eq. (|A2|) can be written as: 



3ff(T - T c ) 
2ir 2 



A 

dfc- 
T 



Ck 2 



r :k 2 



3gT c 

2tt 2 c 



tI(t). (A4) 



The first integral appearing in the previous equation ap- 
proaches a constant as r — > 0, whereas the integral J(r) 
is divergent as r — > as 



/(r) 



-1/2 



(A5) 



To the leading order in r, the left hand side of Eq. 
can be neglected, leading to 



T -T r 



-1/2 



(A6) 



X'^T- 1 ~(T-T C )" 2 . 



(A7) 



This allows to determine the exponent 7: 

7 = 2. (A8) 

Within the self-consistent Hartree approximation, there 
are no k-dependent corrections to the mean field propa- 
gator, all the corrections being included in the renormal- 
ized mass, r. Therefore, 



X (k) = [r + (1 - AnWX 4 )k 2 ] 1 = X [l + (H) 2 ] \ 
where 

^ 2 = (1 - 4itW\ 4 )/t - (T - T c y 2 . 
This immediately leads to 

(All) 



(A9) 
(A10) 



T 1 



Similarly, it is possible to verify that also for W > W c , 
at the spinodal line of the paramagnetic phase, located 
at T = 0, the susceptibility, x(fc m ), and the correlation 
length, £, diverge with the usual Hartree critical expo- 
nents. 



APPENDIX B: FLUCTUATUION-INDUCED 
FIRST ORDER TRANSITION TO THE 
LAMELLAR PHASE 

In this appendix we determine the first order transition 
temperature Tl(W, A) from the homogeneous phase to 
the lamellar phase. We follow the approach taken in (2?| 
for the Coulomb case and we generalize it to the case of a 
finite screening length of the repulsive potential. Let us 
introduce a spatially varying external field, h^, linearly 
coupled to the order parameter field 0k- As a result, 4>k 
is now the sum of an average component, mk = (<ftk), 
and a fluctuation around it £k = 0k — ""ik- The resulting 
equation of state for the average components reads (27| : 



K= r + k 2 



x- 



k 2 



m k 



(Bl) 



+9 



d J ki d J k 2 
(2tt) 3 (2tt) 3 



[mkxTOks + 3G(ki,k 2 )]m k _ kl _ k2 



where the connected correlation function G(k, k') = 
(£ k £ k '), is obtained self-consistently by solving 



TG _1 (k, k') = ro + fc 



A- 2 + /c 5 



5(k + k') (B2) 



3 5 



d 3 q 

(2tt) 3 



m q m k +k'- q + G(q, k + k' - q) 
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together with the unitariety condition: 

/ ^3 G ~ 1 ( k .q) G (^ k ') = ^(k-k'). (B3) 

Note that in the paramagnetic phase, when all /i k , and 
TOk are equal to zero, Eqs. (|B1[) and (|B2p reduce to 
Eq. (HUD with G(k, k') = G(k) S(k + k'). 

In the lamellar phase, characterized by the periodic 
order, instead, we consider: 



and 



h k = h(S(k-k m )+S(k + k m )), (B4) 



m k = m(5(k - k m ) + S(k + k m )) , (B5) 



where k m is given in Eq. (|14p . As shown by Brazovskii 
for a related model [34j |. in this region, the fluctuations of 
wave- vector k with |k| = fc m are dominant, and the effect 
of the off-diagonal terms with k ^ k' can be neglected in 
the correlation function. As a result, 



TG?- X (k) = r + k 2 



A-kW 
A" 2 + k 2 ' 



(B6) 



This expression is formally equivalent to that of Eq. (fTU)) . 
However, in the case of a phase characterized by 
Eqs. (|B4[) and (|B5[) the renormalized mass, r, is given 
by: 



r = r + 3.g 
= r + 3.9 



d J k 

(2tt) 3 
d 3 k 



(G(k) + |m k | 

r 



(B7) 



(2tt)3 r + fc 2 + jij^ 



6g|m| 



By introducing Eqs. (|B4[) - (|B7p in Eq. (|Bip and recalling 
Eq. (fT4"|) . one obtains the following equation of state: 



ft = 



,-2 



r + 2v4ttW - A~ 

., , d 3 k r 

3g 



(2tt)3 r + fe 2 + jfflp 



■3g\m\ 



m 



2V4ttW — A -2 -3ff| 



(B8) 



When we do not consider fluctuations (r — > ro), the last 
equation gives the mean field result. 

Below some temperature, there is a coexistence of the 
paramagnetic phase and the lamellar phase. In zero field 
(h = 0), the former is characterized by m = and the 
latter by m ^ 0, where m is solution of Eq. (|B8[) . i.e., 



(r + 2\/47rIF - A~ 2 - 3.9|m| 2 ) = 0. The transition point, 
which is then associated with a first order transition, is 
obtained as the temperature at which the free energies 
of the two phases are equal. It is convenient to calculate 
directly the free energy difference AF(T) between the 



lamellar (m 7^ 0) and the paramagnetic (m 
at a given temperature T |27| : 



dm 



, OF 
dm 1 



dm'h(m') 



0) phase 



(B9) 



where h(m') is given by Eq. (|B8|) . Changing the inte- 
gration variable from m! to r'(m') solution of Eq. (|B7|) . 
after some algebra we obtain: 



gAF = 



rr(m) I i , 




w*f i_ jU2 _i_ 4-rrtV 



(B10) 



( r l I 1.2 i 4ttW \ ' 

+ ft + A -i +fc i J 



where r(m = 0) is the solution of Eq. (|10p whereas r(m) 
and m are solutions of the following coupled equations: 



r + 2V4^PF - A" 2 - 3g\m\ 2 = 
d 3 k t 



r = r + 3gj^ 



&g\m\ 



(Bll) 



By solving Eqs. (|B10|) and (|B1 1|) numerically for several 
values of W and T (and for g = 1, ro = — 1 and A = 2), 
we have found that, for any given value of W, the sign of 
AF changes at a finite value of T = Tl, that marks the 
first order transition between paramagnetic and lamellar 
phases. The results are shown in Figs. [T] and [H We also 
solved Eqs. ()B9|1 - (|B11|) for several values of W and A 
keeping the temperature T fixed. The transition line in 
the frustration (W / )-range (A) plane is shown in the inset 
of Fig. [T] 



APPENDIX C: SELF-CONSISTENT SOLUTION 
FOR THE DYNAMICS IN THE SPHERICAL 
LIMIT 

In this section we discuss the self-consistent solution of 
the dynamics of the model in the spherical approxima- 
tion. More precisely, our aim is to solve self-consistently 
Eqs. (l2Tj) and (|22|) . In order to do so we remember that: 



Q(t) =r + g 



_cTk_ 

(2tt) 3 



G(k, t). 



(CI) 



From Eq. (|2Tj) . we note that G(k, t) has a maximum at 
a wave vector k m ax(t) given by: 



Q(t) = -2k 2 max (t)t 



4irW\- 2 t 



{^ 2 + k 2 max {t)f 



(C2) 



which gives the expression of Q{t) in terms of the peak 
position k max (t). Inserting Eq. (|C2p into Eq. (|2ip we 
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obtain that G(k,<) = G(k,0) e -M*/(fc,t) ) where 
47rT^fc 2 



/(M) 



A-2 + /c 2 

47rV^A- 2 fc 2 



2fc^^(0 (C3) 



(A- a + *&«(*)) 2 



Obviously, f(k,t) has a minimum for k = k max (t). 

Now we can evaluate the integral in Eq. (|Cip by means 
of the saddle point method, provided that Mt is large. 
This yields: 



Q(t) 



D 



y/Mtf»(k max ,t) 



x exp 



Mt 1- 



4ttVF 



(C4) 
,(*) , 



where B = 2g (7r/2) 3 ^ 2 G(k max , 0). Now, we can relate 
(j(t) to k max (t) through Eq. (|C2[) and, inserting the re- 
sult in the previous equation, we obtain Eq. (|23[) . The 
function P equals: 



P 

where F(k r , 
F(k m 



1 



r 2/c rrtaa . (t) 



4ttVFA- 2 



(A- 



k 2 (+^ 2 



^tkmax rJ^v^max (^)) 



„(*)) 



di 



is defined as: 



m = i - 



4:TtWX~ 2 



(A" 



jL2 

max 



(C5) 



(C6) 



Inserting the right hand side of Eq. (j2"3"|) into the expres- 
sion of the time dependent structure factor we obtain: 



G(k,t) = G(k,0) 



P^MtF{k max {t)) 

k-max 



-MtAf(k,t) 



(C7) 

where Af(k,t) = f(k,t) - f{k max ,t). The position of 
the peak is given by the solution of Eq. (|23|) . Taking the 
logarithm of both sides of Eq. (f23l) and dividing by 
we obtain an asymptotic expression for k max (t) which is 
readily analyzed in the limit t —> oo. 



1 



4ttW 



(*- 2 + kl ax (t)) 2 



•t 4 

ina x 



(*) 



(C8) 



In -r - 2k 2 max (t) ~ 



1 ln[F(fc mtM (t))] 

2 Mt 



hi[fc ma xft)] 
Mt 



where we have neglected the terms which vanish in the 
asymptotic limit, t — > oo. As explained in Sec. IIV1 de- 
pending on the values of W and A, different solutions 
of the previous equation are found. In particular, one 
needs again to distinguish between the cases W < W c 
and W > W c . 



G 



-l 



Go 



-l 




+ 



+ 



FIG. 7: Diagrams of the Self-Consistent Screening Approxi- 
mation Eld, EE El. 



APPENDIX D: SELF-CONSISTENT SCREENING 
APPROXIMATION 

Here we discuss how to compute the complexity within 
the SCSA. This section is unessential, since all the cal- 
culations have been already extensively reported in the 
literature (see Ref. 26]). Nevertheless, we believe that 
this appendix might be useful for the sake of complete- 
ness. 

Following Ref. f§6[ , we first summarize a few properties 
of the matrices in the replica space with structure similar 
to Eq. (j4Tj) . Introducing a matrix E such that E a b = 1 
and the unitary matrix 1 = 6 a b, it is easy to see that the 
product of any two m x m matrices with structure: 

A = ail + a 2 E (Dl) 

is given by: 

A B = (ai&i) 1 + (ai& 2 + a 2 bi + ma 2 b 2 ) E. (D2) 
This relation leads to 



1 



a-2 



a\ ai (ai + 77102) 



■E 



(D3) 



for the inverse of a matrix A. These properties will be 
used in the following. 

In the SCSA, the self-energy, is given by (see 
Fig. [7]): 

S Q6 (k) = ^J ^3 G Qb (k + q)£> o6 (q), (D4) 



where 



D ab (q)= (gT) 1 +n ab (q) 



(D5) 



is determined self-consistently by the polarization func- 
tion 



n afc (q) 



d J p 



G a b(q + p)G fcQ (p). 



(D6) 



The ansatz Eq. (|4ip for the correlation functions implies 
an analogous structure for £ a ;,(k) and n a ;,(k). Inserting 
this ansatz into II a h(q) gives: 



n(q) = [n e (q) - LV(q)] 1 + IV(q)E. 



(D7) 
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where the diagonal and off-diagonal elements of the po- 
larization function are: 



n S (q) = J ^S(q + P )S(p), (D8) 
IV(q) = J ^^(q + p)^(p). 



Using the property expressed in Eq. (|D2|) . it is now 
straightforward to determine -D a &(q) which, in the limit 
to — > 1, is given by: 



Dg(q) = [Dg(q) - 2V(q)] 1 + 2V(q)E, 



where 



Vg(q) = [(gT) 1 +n g (q) 
IV(q)2?2 (q) 



(D9) 



(D10) 



l-IV(q)2? g (q) 



Analogously inserting the above equation into Eq. (|D4jl . 
we get for the self-energies 



E(q) = [E„(q) - E^(q)] 1 + E^(q)E, 



(Dll) 



where 



E e (q) 
S^(q) 



d 3 p 

TV 7 (2tt) 3 
2 /" d 3 p 



TV 7 (2tt) 3 



2>s(p)S(q + p), (D12) 
P^(p)^(q + p). (D13) 



This set of equations is closed by the Dyson Equation, 
Eq. (gOl): 

G- b \k) = [G (k) + E s (k) - E^(k)] S ab + E^(k), 

(D14) 

which, in the limit m — > 1, according to Eq. (|E)3|) . gives: 



^- 1 (k) = G - 1 (k) + E e (k) 
for the diagonal elements, and 

5 2 (k)E^(k) 



= S(K - 



1 - 0(k)E,r(k) 
1 



(D15) 



(D16) 



£-i(k)-E^-(k) 



for the off-diagonal ones. 

Within the self-consistent screening approximation, 
the free energy assumes the following expression [l8|, [l9l . 
[11 HI: 



F(m) 
2toT 



= TrlnG" 



TrlnZT 



TrGE. 



(D17) 



The evaluation of the traces in the replica space is 
straightforward. After that we are able to perform the 
derivative with respect to the number of replicas in order 
to compute the complexity, E, according to Eq. (l38|) . To 
this aim we have to take the analytical continuation to 
to — > 1 . The configurational entropy can be written as 
the sum of two contributions [2a, um : 



E = E« 



,(2) 



(D18) 



where 



d J k 

(2tt)3 



In 1 - 



C7(k); + C7(k)/' 



(D19) 



and 



In 1 



d J k 

(2tt) 3 

grn^(k) 
i + . g rn e (k) 



grn^(k) \ 
i + 5rn e (k)y 



(D20) 



From Eqs. (|D19IID20p and Eq. (|Dl6j) . it immediately fol- 
lows that E = if ^(k) vanishes or, equivalently, if 
Ejr(k) vanishes. 



Solving numerically Eqs. (|D8|) - (|D16[) we get the ex- 
pressions of the correlators, £/(k) and ^"(k), of the self- 
energies, Eg(k) and E^-(k), and of the polarization func- 
tions, Ilg(k) and Iljr(k). Then, according to Eqs. (|E)19|) 
and (1D20[) we are able to compute the complexity, E. 
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